function prior = conjugateprior(D, modelname, place_prior_on_mixing)


if exist('place_prior_on_mixing', 'var') && place_prior_on_mixing
    prior.a0 = 0.001;
end

prior.m0 = zeros(1, D);

switch modelname
    case 'gmm'
        prior.b0 = 0.001;
        prior.W0 = eye(D)*1000;
        prior.v0 = D;

    case 'kmeans'
        prior.T0 = eye(D)*1000;
 
    otherwise
        error('no model selected');
end
        
